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Abstract. We present an overview of our recent work on transport phenomena 
in Weyl metals, which may be connected to their nontrivial topological properties, 
particularly to chiral anomaly. We argue that there are two basic phenomena, which 
are related to chiral anomaly in Weyl metals: Anomalous Hall Effect (AHE) and 
Chiral Magnetic Effect (CME). While AHE is in principle present in any ferromagnetic 
metal, we demonstrate that a magnetic Weyl metal is distinguished from an ordinary 
ferromagnetic metal by the absence of the extrinsic and the Fermi surface part of the 
intrinsic contributions to the AHE, as long as the Fermi energy is sufficiently close to the 
Weyl nodes. The AHE in a Weyl metal is thus shown to be a purely intrinsic, universal 
property, fully determined by the location of the Weyl nodes in the first Brillouin zone. 
In other words, a ferromagnetic Weyl metal may be thought of as the only example 
of a ferromagnetic metal with a purely intrinsic AHE. We further develop a fully 
microscopic theory of diffusive magnetotransport in Weyl metals. We derive coupled 
diffusion equations for the total and axial (i.e. node-antisymmetric) charge densities 
and show that chiral anomaly manifests as a magnetic-field-induced coupling between 
them. We demonstrate that an experimentally-observable consequence of CME in 
magnetotransport in Weyl metals is a quadratic negative magnetoresistance, which 
will dominate all other contributions to magnetoresistance under certain conditions 
and may be regarded as a smoking-gun transport characteristic, unique to Weyl metals. 
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The last decade has seen an explosion of research in topologically-nontrivial states 
of matter, following the remarkable discovery of topological insulators (TI) [TJ [21 
El IH El El 0 i- One normally associates topologically nontrivial properties with 
insulators, the spectral gap providing the rigidity and insensitivity to fluctuations, 
which are characteristic of these states of matter. One of the most recent notable 
developments in this field, however, has been the realization that even gapless metallic 
and semimetallic states may be topologically nontrivial in much the same sense as 
gapped insulators mmmm- These recent developments have partly been anticipated 
in the earlier pioneering work of Volovik [[13], promoting topological classification of all 
possible fermionic ground states, which shed new light on such wide-ranging phenomena 
as the robustness of the Fermi liquid and the hierarchy problem in particle physics. 

The recent work focussed on Weyl semimetals [a nm nn H2], which, in Volovik’s 
classification, belong to the Fermi point universality class of fermionic vacua. In 
the condensed matter context a Weyl point or node is a point of contact between 
two nondegenerate bands in the first Brillouin zone (BZ). Such electronic structure 
features have in fact been noticed and studied since the earliest days of the theory of 
solids mmm, but only recently have their topological properties come into focus. 
Two important ingredients in the appearance of Weyl node features in the electronic 
structure are broken time reversal (TR) or inversion (I) symmetries and dimensionality 
of space. The broken symmetry requirement is a consequence of the well-known Kramers 
theorem. If both TR and I symmetries are present, the band eigenstates in a solid must 
be (at least) doubly degenerate at every value of the crystal momentum in the first 
BZ. Contact can then occur only between pairs of bands, which is in general impossible 
without fine-tuning. The contact between two nondegenerate bands is, however, possible 
generically in three spatial dimensions (3D). Indeed, near the point of contact, the 
momentum-space Hamiltonian must have the following form, dictated by the so-called 
Atiyah-Bott-Shapiro construction HU 

H = ±uf<t • k, (1) 

where the triplet of Pauli matrices a describes the pair of touching bands, the sign in 
front refers to two possible chiralities of the band contact point, and we have subsumed 
any possible spatial anisotropies into the definition of the crystal momentum k. We will 
use h — c — 1 units throughout this paper, except in some of the final results. The 
point of band degeneracy occurs when all three components of the crystal momentum 
vanish. This is why three-dimensionality is important: we must have three momentum 
components available as “tuning parameters” to create a point of degeneracy between 
the two bands. A naively analogous degeneracy point in 2D graphene in fact does not 
exist, as was pointed out in the seminal work of Kane and Mele [T|, which started the 
field of TI. It only looks like a degeneracy point due the smallness of the spin-orbit (SO) 
interactions in graphene, which always open up a small gap and destroy the degeneracy. 
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The presence of Weyl degeneracy nodes is, in principle, a common feature of the 
electronic structure of all 3D magnetic or noncentrosymmetric materials. In most cases, 
however, their presence is of no consequence, since they will generically occur far away 
from the Fermi energy [18]. Weyl semimetal refers to a situation in which the Fermi 
level exactly coincides with the Weyl nodes and no other states are present at the Fermi 
energy. This situation is not generic and occurs only under special, but not entirely 
uncommon circumstances. 

To identify these, it is convenient to start from the situation in which both TR 
and I symmetries are present. TR symmetry requires that if a Weyl node is present at 
momentum k, another node with the same chirality must be present at — k. I symmetry, 
on the other hand, requires that if a Weyl node is present at k, a node of opposite 
chirality must be present at —k. This implies that in the presence of both TR and I, 
there may exist a pair of opposite-chirality Weyl nodes at the same TR and I-invariant 
crystal momentum. For concreteness, we will take this momentum to be at the BZ 
centre, i.e. at the T point. The minimal k • p Hamiltonian, describing this situation, 
has the form 

H = v f t z (t ■ k + A r x , (2) 

where r is an additional set of Pauli matrices that describes the two nodes of opposite 
chiralities, i.e. H has the form of the Hamiltonian of a 3D four-component Dirac fermion 
with “mass” A. The mass term A annihilates the two Weyl nodes, unless it vanishes 
by symmetry or is fine-tuned to zero. To realize the former situation, the states at 
the T point must either form a four-dimensional irreducible representation of the space 
group of the crystal, as proposed to occur in /3-cristobalite Bi0 2 [IS] and, with a few 
modifications (Dirac Hamiltonian is replaced by Luttinger Hamiltonian in this case), 
in some of the pyrochlore iridate materials (2HI Si 122] - Or they must belong to two 
distinct two-dimensional representations, so that pairwise crossing, protected by crystal 
symmetry, is possible, as realized in the recently discovered Dirac semimetal materials 
Na 3 Bi and Cd 2 As 3 [23] E3 [25l EH EZj . 

The latter situation, i.e. with A fine-tuned to zero, occurs at the transition point 
between an inversion-symmetric 3D TI and a normal insulator (NI) [251 fill 1251 ISO] . In 
the present review we will focus exclusively on this case. Our results, however, do not 
depend on this choice and are applicable to all realizations of Dirac and Weyl semimetals. 
As mentioned above, Weyl semimetal may be obtained from Dirac semimetal by breaking 
either TR or I. We will focus on the TR-breaking case here, since it produces the simplest 
possible kind of Weyl semimetal, with only a single pair of Weyl nodes. Generalization 
to many pairs is in most cases straightforward, but it does complicate the technical 
details of course. Formally, this type of Weyl semimetal is obtained by adding a Zeeman 
spin-splitting term to the Dirac Hamiltonian Eq. [2] 

H = v F r z cr ■ k + At x + bcr z , (3) 

where the spin-splitting term may arise either from an external magnetic field, or 
from an intrinsic magnetic ordering, and we have chosen the ^-direction to be the 
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magnetization direction. To analyze this Hamiltonian, it is convenient to perform the 
following canonical (i.e. commutation relation preseving) transformation of the cr and 
r operators 

a ± —> r z a ± , t 1 * 1 —> a z r ± . ( 4 ) 

The Hamiltonian Eq. [3] then transforms to 

H = vf{c r x k x + cr y ky) + (b + vpk z r z + Ar x )a z , (5) 

Diagonalizing the vpk z r z + At x block of the Hamiltonian we obtain 

H± = vp(cr x k x + a y k y ) + m±(k z )a z , (6) 

where 

m±{k z ) =b± \Jv 2 F k 2 z + A 2 . (7) 

Eq. [b] looks like the Hamiltonian of a pair of 2D Dirac fermions with “masses” m±(k z ), 
which depend on the ^-component of the crystal momentum as a parameter. The mass 
m + (k z ) is always positive, but m-{k z ) may vanish and change sign if b > A. If this is 
the case, vanishes at k± = ±k 0 , where 


1 


ko = — Vb 2 - A 2 . 


vf 


( 8 ) 


These points along the £-axis in momentum space, where the “mass” m_(/c 2 ) vanishes, 
are the locations of the Weyl nodes. The two bands, which are the eigenstates of the 
block of the Hamiltonian, are degenerate at those points. Since massive 2D Dirac 
fermions are associated with half-quantized Hall conductivity 

„2 


a 


2D 


= —sign (m), 


(9) 


( 10 ) 


*» 2ft‘ 

it follows that the Weyl semimetal at b > A has a Hall conductivity given by 
e 2 

&xy "^“2/c 0) 

i.e. is equal to the distance between the Weyl nodes in units of e 2 /h [22]. 

It is easy to show that this Hall conductivity may be associated with chiral edge 
states. Indeed, suppose the sample is finite in the ^-direction and we will take the 
sample to occupy the y < 0 half-space. Replacing k y —* —id/dy, the H block of the 
Hamiltonian becomes 


d 

H_ = — ivp-p- & v + Vpcr x k x + rri-(k z , y)cr z ■ 
dy 


( 11 ) 


Let us first set k x = 0 and look for a zero-energy solution of the Schrodinger equation 
(SE) 

H-'&{k z ,y)= 0, (12) 

localized at the sample boundary y — 0. We look for a solution in the form 
'ft(k z ,y) = ia y e^ kz,yS> 4> : 


( 13 ) 
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where f(k z ,y) is a scalar function, while 0 is a two-component spinor. Plugging this 
into the SE Eq. 12, we obtain 
df{k z ,y) 


v F - 


dy 


+ m-(k z ,y)a x 


= 0 . 


(14) 


The solution of Eq. [14} satisfying our requirements is 

1 


f(kz,y) = — m-{k z ,y')dy', 


(15) 


vf Jo 

and 0 = \a x = —1), i.e. the eigenstate of <j x with eigenvalue —1. Thus we finally obtain 
the following result for the solution of Eq. [14} localized at the sample boundary 


*(k z ,y) = = 1). (16) 

This solution exists as a localized edge state as long as m_(k z ,y —)• — oo) > 0, i.e. as 
long as — ko <k z < ko . For this reason it is called a Fermi arc [9j. It is straightforward 
to see that the dispersion of this state in the ^-direction is given by e = v F k x , i.e. it is 
chiral. 

At this point we will wrap up the introductory part of this review and move 
on to a more detailed account of the transport properties of Weyl semimetals and 
metals, based on a more realistic microscopic model of Weyl semimetal. The 
rest of the paper is organized as follows. In Section [2] we introduce a realistic 
microscopic model of a Weyl semimetal, based on a magnetically-doped TI-NI multilayer 
heterostructure. In Section [3] we discuss basic electromagnetic response properties of 
this model Weyl semimetal and introduce two distinct components of the response, 
which are directly related to the nontrivial topology of Weyl nodes: the Anomalous 
Hall Effect (AHE) [33;, [3TJ and the Chiral Magnetic Effect (CME). In Section [4] we 
extend the theory of electromagnetic response of Weyl metals to diffusive transport 
regime. We demonstrate that magnetic Weyl metals are distinguished from ordinary 
ferromagnetic metals by a lack of impurity-scattering and Fer mi -surface contributions 
to their anomalous Hall conductivity. Instead their anomalous Hall conductivity is 
determined by the distance between the Weyl nodes in momentum space and nothing 
else. We also show that CME manifests in the diffusive regime as a new kind of 
weak-field magnetoresistance, which is negative and quadratic in the field. This novel 
magnetoresistance is expected to occur in all types of Weyl metals and may be regarded 
as their smoking-gun transport characteristic. We conclude in Section [5] with a brief 
overview of the main results. 


2. Heterostructure model of Weyl semimetal 

As discussed in the Introduction, perhaps the most straightforward (at least 
theoretically) way to realize a Weyl semimetal phase is to break TR symmetry in a 
material, which in its nonmagnetic state is naturally poised near the phase transition 
between a TI and NI. One way to achieve this is to engineer a composite material, made 
of alternating thin layers of TI and NI, as shown in Fig. [lj This system may be regarded 







Chiral anomaly and transport in Weyl metals 


6 



Figure 1. Cartoon of the heterostructure model of Weyl semimetal. Magnetized 
impurities are shown by arrows. 

as a “hydrogen atom” of Weyl semimetals: the most elementary yet realistic system, the 
description of which is simple enough that purely analytical theory of many phenomena 
is possible, as will be seen below. 

This system may be described as a chain of 2D Dirac surface states of the TI 
layers, which are coupled by tunnelling matrix elements A s for a pair of surface states, 
belonging to the same Tl layer, and Ap for a pair belonging to nearest-neighbor Tl 
layers. We will assume, for concreteness, that A $,d A 0. The corresponding Hamiltonian 
has the form 

H = v f t z (z x cr) ■ k ± 8 itj + A s r x 5 id + ^ A D (r + 8 jti+1 + t - <^_i) , (17) 

where i,j label the unit cells of the superlattice in the growth (A direction, ky = (k x , k y ) 
are the crystal momentum components, transverse to the growth direction, Pauli 
matrices r act on the which surface pseudospin degree of freedom, while cr act in the spin 
degree of freedom. Diagonalizing the growth direction hopping part of the Hamiltonian 
by Fourier transform and performing the canonical transformation of Eq. [4j we obtain 

H = v F (z x cr) • k + A (k z )a z , (18) 

where 

A (k z ) = A s t x + (r + e ikzd + h.c.) , (19) 

d being the period of the superlattice in the growth direction. A (k z ) may now be 
diagonalized separately, as it commutes with the Hamiltonian. Its eigenvalues are given 
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A (k z ) 



+ A 2 D + 2A s A d cos (k z d). 


( 20 ) 


The corresponding two-component spinor wavefnnctions are 

/ A _ i A _ „—ikzd ' 


\u (k z )) = 7= 


V2 V 


(i,t As + Ace " 




( 21 ) 


where t = ± labels the two eigenvalues. The spin block of the Hamiltonian may now be 
diagonalized as well, giving rise to the following eigenvalues 


e s t(k) = se t (k) = s^v 2 F {k 2 x + k 2 ) + m 2 t {k z ), 
where s — ± and 


( 22 ) 


m t (k z ) = tA(k x ). 


(23) 


The double degeneracy of the band eigenvalues at every k in Eq. [22] is the Kramers 
degeneracy, arising due to the presence of both TR and I symmetries. The corresponding 
eigenvectors are given by 


,st 


( k )) _ \ \l 1 + s 


m t {k z ) 

et(k) 


-ise 



m t (k z 
' et(k) 


(24) 


where e 1 ^ = -^^==. The full four-component eigenvectors of the Hamiltonian Eq. 
may be represented as a tensor product ^(kg)) and |n st (k)): 
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^(k)) = |^(k))0|^(fc 2 )). 


(25) 


The heterostructure, described by Eq. [17] can exist in two distinct phases: strong 
3D TI and an NI. This can be easily seen by computing the index, using the method 
of Fu and Kane |35j . Namely, we compute the eigenvalues of the parity operator t x (or 
T x a z after the canonical transformation of Eq. [4]) at two TR-invariant momenta in the 
first BZ: Ti = (0, 0, 0) and T 2 = (0, 0, n/d). Taking the product of the parity eigenvalues 
over all the occupied bands and over r 12 , the Z 2 index is found to be 

(-1)" = sign(A 5 - A d ). (26) 


Thus the multilayer is a strong 3D TI when Ap > A$ and an NI otherwise. The point 
A s = Ajj marks the transition point between the 3D TI and NI. At this point, the gap 
at r 2 closes and the heterostructure is a Dirac semimetal. 

To obtain a Weyl semimetal, we break TR symmetry in the heterostructure, 
assuming that it is tuned close to the TI-NI transition point and the band gap |Ag — Ad| 
is not large. The TR breaking may be accomplished by doping the heterostructure 
with magnetic impurities. At sufficient concentration, the impurities will form a 
ferromagnetic (FM) state, and we will assume that the magnetization points along the 
growth direction of the heterostructure (distinct, non-Weyl nodal states are obtained 
if the magnetization is in the xy-plane [36]). We will describe the FM state of the 
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heterostructure by adding a term ba z 5ij to the Hamiltonian Eq. 17 
degeneracy of the t — ± Kramers doublet, modifying Eq. 23 as 

m t (k z ) = b + tA(k z ). 


This splits the 


(27) 


Taking the spin splitting b to be always positive, the “mass” m + (k z ) is then always 
positive, while m_{k z ) may vanish and change sign as a function of k z , if the spin 
splitting is sufficiently large. This happens when b > b c \ = A 5 — Ajj\- The Weyl node 
locations along the k x = k y = 0 line in the crystal momentum space are given by the 
solutions of the equation 


m~(k z ) = b - A(k z ) = 0, 
which gives kf = Ti/d ± k 0 , where 


ko = - arccos 
d 


A| + A 


f 


(28) 


(29) 


2A s A d 

When the spin splitting reaches the upper critical value b c 2 = As + Ad, ko = n/d 
and the Weyl nodes are annihilated again at the BZ edge. The resulting state is a 3D 
quantum anomalous Hall insulator with a quantized Hall conductivity 

Q 2 


@xy 


e 

hd 


(30) 


3. Electromagnetic response of Weyl semimetals 


3.1. General remarks 


Topologically-nontrivial states of matter, such as 3D TI, and Weyl semimetals, have a 
distinguishing spectroscopic feature: the presence of metallic edge states. These surface 
metallic states are unusual and may only exist on the surface of a bulk 3D topological 
phase. As discussed in the Introduction, in the Weyl semimetal state, of interest to us 
here, the Fermi surface of this metal is not a closed 2D curve, as it must be in any regular 
2D metal, but instead forms an open Fermi arc, with the end points at the locations of 
the bulk Weyl nodes, projected onto the surface BZ plj f351J . 

However, apart from such spectroscopic distinguishing features, topological phases 
often have unusual electromagnetic response characteristics, which are always of 
particular interest, as these are the manifestations of the unique quantum physics of 
these phases on macroscopic scales. In 3D TI, this unusual electromagnetic response 
may be described by the so-called 0 -term in the action of the electromagnetic field [SJ 


S e = 


e 2 6 
32 F 2 


e 2 6 !' 

dt d 3 re> wap F llv F a a = — / dt d 3 rE ■ B. 


47T 2 


(31) 


The parameter 6 is equal to 7 r in 3D TI, the only nonzero value, compatible with TR 
symmetry jHj. The 0-term in Eq. 31, however, is a full derivative, and thus has no 
effect on Maxwell’s equations in the bulk of the TI. Its only real effect is to generate 
the half-quantized AHE on the sample surfaces in the presence of surface magnetization 
(e.g. due to magnetic impurities). 
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When TR is broken in the bulk, as it is in the Weyl semimetal state, a 0-terrn with 
non-constant 6 is now allowed by symmetry 

Se = ^J dt d 3 r0(r)e^“ fl F /w F Q3 . (32) 

The only functional form, compatible with translational symmetry, is linear, i.e. 

6( r) = 2b ■ r, (33) 

where b is a vector, which is odd under TR and even under 1 and coincides, as will be 
seen explicitly below, with the bulk spin splitting. If, in addition to TR, I symmetry is 
violated as well, the “axion field” 6 may also acquire a linear time dependence 

6{ r, t) = 2b • r — 2 b 0 t, (34) 


where bo is the energy difference between the Weyl nodes, allowed if / is violated. Unlike 
the 0-term in 3D TI, Eq. 32 does modify Maxwell’s equations in the bulk of the Weyl 
semimetal and thus has observable consequences, namely the already mentioned AHE 
and the Chiral Magnetic Effect (CME) [39], both of which will be discussed in detail 
below. 

The action of Eq. [32] may be described as being a consequence of chiral anomaly, 
an important concept in relativistic field theory [40, j4I| 02], which has recently found 
its way into condensed matter physics and has been realized to play an important role 
in the theory of topologically-nontrivial states of matter [43, 44] 03J. To gain a basic 
understanding of chiral anomaly, it is useful to go back to the generic low-energy model 
of Weyl semimetal, given by Eq. [3} For simplicity, we will subsume the Fermi velocity 
Vf in Eq. [3] into the definition of momentum and set the “mass term” A = 0. Then we 
obtain 


H = T z cr ■ k + b 0 r z + b • <r, 


(35) 


where we have included an energy difference between the Weyl nodes, described by the 
b 0 r z term, which is allowed when the I symmetry is violated, and oriented the TR- 
breaking vector b along an arbitrary direction. It is convenient to represent the system 
in terms of an imaginary time action, including possible coupling of the electrons to an 
electromagnetic field 

S = f drd 3 r [d T + ieA 0 + b 0 r z + r z cr ■ (—IV + eA + hr z )] ip, (36) 


where A M = (A 0 ,A) is the electromagnetic gauge potential and ip\ip are the 4- 
component spinor Grassman fields. We have suppressed all explicit spinor indices in the 
Grassmann fields for brevity. We now observe that in addition to the charge conservation 
symmetry, the imaginary time action Eq. 36 possesses an extra chiral symmetry 

ip e- irZ0/ V, (37) 


which expresses an apparent separate conservation of the number of fermions of left and 
right chirality. This suggests that the terms boT z and b • cr in Eq. [36] could be eliminated 
by a gauge transformation 

Ip -> ^ ^ e ir‘tl( r ,T)/2' 


( 38 ) 
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where 9( r, r) = 2b • r — 2ib$r and one should keep in mind that 0 and -01 are not 
complex conjugates of each other, but are independent variables in the fermion path 
integral. The imaginary time action then becomes 


S= drdr {d T + ieA 0 + t z ct ■ (—iV + eA)] -0, 


(39) 


which describes two Weyl nodes of opposite chirality, existing at the same point in 
momentum space and in energy. This argument then leads one to the conclusion that 
the system of Weyl nodes, separated in energy and momentum, is equivalent to the 
system of two degenerate Weyl nodes and thus does not possess any special transport 
properties, which is incorrect. What is missing in the above argument is exactly the 
chiral anomaly: while the imaginary time action Eq. [36] does indeed possess the chiral 


symmetry, the gauge transformation of Eq. 38 changes not only the action itself, but 
also the measure of the path integral. This change of the path integral measure is what 
gives rise to the 0-term Eq. [32] This may be shown explicitly using the Fujikawa’s 


method uni STj. We refer the reader to Ref. m for technical details. 

To see the physical consequences of the 0-term, we integrate Eq. 32 by parts and 
eliminate a total derivative term. This leads to the following action, which resembles 
the Chern-Simons term in 2 + 1-dimensions 

S e = / dtd 3 r d,9e^A v d a Ap. (40) 


Varying Eq. 40 with respect to the vector potential, we obtain the following expression 
for the current 


3, = “dcAf, 


A = 1,2,3, 


and 


Jv = 


2tt 2 


boe 0ual3 d n A 


aSi/3- 


(41) 


(42) 


Eq. 41 clearly represents the AHE, while Eq. 42 represents another effect, related to 
chiral anomaly in Weyl semimetals, the Chiral Magnetic Effect [39], whose physical 
meaning is somewhat subtle and will be discussed in detail later. 

However, the picture of chiral anomaly and related transport phenomena, presented 
above, while simple and appealing, is not fully satisfactory, for the following reasons. 
Chiral anomaly is a sharply-defined concept in the context of relativistic held theory, 
where massless fermions in unbounded momentum space possess exact chiral symmetry 
Eq. 37, which may be violated by the anomaly. In the condensed matter context, 
however, the situation is less clear. Even though chiral symmetry may be formally 


defined in a low-energy model of a Weyl semimetal Eq. 35, in which the band dispersion 
is approximated as being exactly linear and unbounded, no real microscopic model of 
Weyl semimetal actually possesses such an exact symmetry. It may only appear as an 
approximate low-energy symmetry. Since the chiral symmetry is not present to begin 
with, it is then unclear to what extent is it meaningful to speak of its violation by 
chiral anomaly. The role of impurity scattering, present in any real condensed matter 
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system, and important even conceptually in any discussion of transport phenomena; 
as well as the role of finite electron or hole density, present when the Fermi energy 
does not exactly coincide with the location of the Weyl nodes, are also completely 
unclear in the above field-theory discussion. In the remainder of this paper we will thus 
develop a fully microscopic theory of chiral anomaly and its manifestations, including the 
contribution of impurity scattering, and finite charge carrier density, based on the TI-NI 
heterostructure model, described in Section [2j which does not possess chiral symmetry, 
yet exhibits both the AHE and CME, described above. 


3.2. Microscopic theory: Anomalous Hall Effect 


We will now go back to our microscopic model of a Weyl semimetal, described in 
Section [2} We want to explicitly derive the topological term in the action of the 
electromagnetic held by integrating out electron variables, coupled to the held. We 
then start from the imaginary time action of electrons in the Weyl semimetal, coupled 
to electromagnetic held 

S — J drd 3 r j'h^r, r) d T — p + ieA 0 (r,r) + H T(r,T)j, (43) 

where A 0 (r, r) is the scalar potential and 

H = vpT z (z x <j) ■ (— i'V + eA) + A + ba z , (44) 

is the Hamiltonian of noninteracting electrons in Weyl semimetal, minimally coupled to 
the vector potential A. We will ignore the ^-component of the vector potential as it will 
not play any role in what follows. 

Using the results of Section [2] for the eigenvalues and eigenvectors of the Weyl 
semimetal Hamiltonian in the absence of the electromagnetic held, we can now integrate 
out electron variables in Eq. [43] and obtain an effective action for the electromagnetic 
held, induced by coupling to the electrons. This action will contain two distinct kinds 
of contributions. The hrst kind will contain terms, proportional to E 2 and B 2 , where E 
and B are the electric and magnetic helds. These terms describe the ordinary electric 
and magnetic polarizability of the material, and we will not discuss this part of the 
response here. The second kind contains the “topological” contribution, which has the 


form of the “3D Chern-Simons” term of Eq. 40 Adopting the Coulomb gauge V-A = 0, 
we obtain 

5 = J2 e z0a/3 n(q, ift) A 0 (-q, -in)q a A p { q, ifi), (45) 

where q a = q a /q and summation over repeated indices is implicit. The ^-direction in 


Eq. 45 is picked out by the magnetization b. The response function n(q, iVt) is given by 

%[6't'(k)] - n F [i s t {k + q)] 


n(q, fD) = 


ie 2 VF 


V 


E 


x 


(~k+qUl 


k 

s’t' 


iVt + £ s 't'(k) — £ s t(k + q) 

>(*k*V Wl4+ q >, 


(46) 
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where £ st (k) = e s t(k) — and summation over repeated band indices s,t is implicit. 

To evaluate II (q, iQ) explicitly it is convenient to rotate coordinate axes so that 
q = qx and assume that ep > 0. The ep < 0 result is evaluated analogously. As seen 
from Eq. 46, there exist two kinds of contributions to the response function II (q, ifl): 
interband, with s ^ s', and intraband, with s = s' = +. Both in general contribute and 
need to be taken into account. Let us first evaluate the interband contributions. In this 


case we can set iCt = 0 in the denominator of Eq. 46 from the start. Expanding to first 
order in q and integrating over the transverse momentum components k XiV , we obtain 

rn/d 

/ dk z sign [m t {k z )\ 

' —iv/d 


n 


inter , 


<^>=Se 

t 

\m t (k z )\ 


x 


1 - 


1 - 




0(e F - \m t (k z )\) 


(47) 


where @(x) is the Heaviside step function. As clear from Eq. 47, n mter (q ,ifl) consists of 
two distinct terms. The first term inside the curly brackets in Eq. [47] is the contribution 
of the completely filled s = — bands. The second term is the contribution of incompletely 
filled s = + bands. 


To understand the meaning of Eq. [47} let us now evaluate the intraband contribution 

It is then clear that, unlike 


to n(q, In this case we have s = s' = + in Eq. 46 


in the case of the interband contribution, evaluated above, the value of the intraband 
contribution depends on the order in which the limits of q —)■ 0 and O —y 0 are taken. 
If the limit is taken so that fi/up|q| —>■ oo, then one is evaluating the DC limit of 
a transport quantity, i.e the optical Hall conductivity. In this case, the intraband 
contribution vanishes identically. On the other hand, if the limit is taken so that 
D/V|q| —> 0, then one is evaluating an equilibrium thermodynamic property, whose 
physical meaning will become clear below. In this case, the intraband contribution is 
not zero and is given by 

dnp(x) 


n 


intra 


(q, «n) = 


ie 2 VF 

V 


EE 


dx 


(4+ql-k X*k Kel^k+q)- 


x=et(k )—ep 


(48) 


The derivative of the Fermi distribution function in Eq. |48| expresses the important fact 
that H intra (q, ifi) is associated entirely with the Fermi surface, unlike n mter (q, *f2), to 
which all filled states contribute, including states on the Fermi surface. Expanding to 
first order in q as before and evaluating the integrals over k x . y , we obtain 


n mtra (q,iD) = 


e 2 q 
8i r 2 


E 


t /d 


dk- 


m t (k z 


-0(e F - \m t (k z )\). 


(49) 


l ir/d 

i.e. the intraband contribution to n(q, iit) is equal to the second term in the square 
brackets in Eq. 47 in magnitude, but opposite in sign. Combining the inter- and 
intraband contributions to n(q,if2) we thus obtain 

t /d 


P i { 7 w _ r n / a 

= 8^E f • 


dk z sign [m t (k z )] [1 - 0(e F - \m t (k z )\)\, (50) 
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(a) 


cr xy h/e 2 



Figure 2. (a) Plot of the band edges along the z-direction in momentum space. The 
parameters are such that two Weyl nodes are present, (b) Total intrinsic anomalous 
Hall conductivity (solid line), a^ y (dashed line), and a^ y (dotted line). Note that the 
van Hove-like singularities in a xy and a xy . associated with band edges, mutually cancel 
and the total Hall conductivity a xy is a smooth function of the Fermi energy. 



o- xy h/e 2 



e F/^S 


Figure 3. (a) Plot of the band edges along the z-direction in momentum space. The 
spin splitting is not large enough for the Weyl nodes to appear (i.e. b <b c i) and the 
spectrum has a full gap. (b) Total intrinsic anomalous Hall conductivity (solid line), 
< 7 xy (dashed line), and a^ y (dotted line). 


i.e. the last term in Eq. 47 cancels out when the low-frequency, long-wavelength limit 
is taken in such a way that h2/fi?|q| —> 0. On the other hand, when h2/ui?|q| —> oo, the 
intraband contribution vanishes and 


n(q,ift) = n mter (q,fO). (51) 

This physical difference in the kind of response the system exhibits is the basis of Streda’s 
separation of contributions to the Hall conductivity into a xy and a r J y [48j. Our analysis 
makes it clear that this is the most physically-meaningful separation of contributions to 
the intrinsic anomalous Hall conductivity, since it corresponds to distinct and, at least 
in principle, separately measurable, contributions to the response function n(q, ifi). 
Generalizing the above results to arbitrary sign of we finally obtain 

S = -ie* 0a/3 a xy f d 3 rdTA 0 (r, r)d a A e (r. r), 


(52) 
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where, if the low-frequency long-wavelength limit is taken so that 0/np|q| —> oo: 

e 2 v-^ [ n ^ d 

cr xy = / dk z {sign[m t (Ax)] [©(e^ + \m t {k z )\) - ©(e^ - \m t {k z )\)] 

O' 71 '*' t J-Tt/d 

+ m ^^ &(\e F \ ~ \m t (k z )\)X . (53) 

\ e F\ J 

This expression corresponds to the full DC anomalous Hall conductivity. On the other 
hand, when the low-frequency, long-wavelength limit is taken so that Q/vp|q| —» 0, we 
obtain 

e 2 v-^ f^ d 

Gxy = a xy= J /d dkz si § n K(^)] 

x [0 (cf + \m t (k z )\) - 0(e F - \m t (k z )\)} . (54) 

This is precisely the Streda’s a xy contribution to the Hall conductivity, which is a 
thermodynamic equilibrium quantity, equal to 


a Z = e 


dN \ 
dB J 


(55) 


where N is the total electron number. This relation follows immediately from Eq. 52 


and the order of limits D/up|q| —y 0, which corresponds to thermodynamic equilibrium. 
Correspondingly, the a T xy contribution is given by 


a xy — a xy 


U _ 
°xy 


8tt 2 ^ 
t 


—ir/d 


\ e F\ 


-@(|e^| - \m t (k z )\). 


(56) 


As clear from the above analysis, cr 1 is the contribution to a xy that can be associated 
with states on the Fermi surface. This contribution is nonuniversal, i.e. it depends on 
details of the electronic structure and is a continuous function of the Fermi energy. a xy , 
on the other hand, is the contribution of all states below the Fermi energy and is a 
thermodynamic equilibrium property of the ferromagnet. It attains a universal value, 
which depends only on the distance between the Weyl nodes, when the Fermi energy 
coincides with the nodes, i.e. when €f — 0 

n ii = tl Q - 

°xy 4?r 2 ’ 

where 



K, — - arccos 
d 


A 2 , + A l-ff 


(58) 


2A5A/) 

is the distance between the Weyl nodes. When b > b c 2 , the Weyl nodes annihilate at the 
edges of the Brillouin zone and a gap opens up. In this case /C = 2tt/ d, i.e. a reciprocal 
lattice vector and a r J y is quantized as long as the Fermi level is in the gap [37] . Both 
contributions, along with the total anomalous Hall conductivity a xy are plotted as a 
function of the Fermi energy in Figs. il in two different cases: when Weyl nodes are 
present and when they are not. The former occurs when b c \ < b < b c 2 . As can be seen 
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from Fig. [ 2 J Weyl nodes provide the dominant contribution to a T J y and to the total Hall 
conductivity a xy , if the Fermi level is not too far from the nodes. 

It is interesting to note the following property, which is evident from Fig. [2j Both 
the total anomalous Hall conductivity a xy and the two distinct contributions to it, cr xy r 
appear to exhibit a quasi-plateau behavior when e F is not too far from the Weyl nodes. 
To understand the origin of this behavior, consider the derivative of a xy with respect to 
the Fermi energy 

da 11 e 2 f 7r / rf 

/ = / dk z sign[m t (k z )\8(e F - \m t {k z )). (59) 

de F 87t J_ n/d 

This is straightforward to evaluate analytically and we obtain 
do 11 e 2 r /d 

= - — / dk z [S(A(k z ) - b + e F ) - <5(A(k z ) - b- e F )] 
oe F 87 t 2 J_n/ d 

= ( 1 /vf+ ~ 1 /v F ~) , ( 60 ) 

where we have assumed that e F is sufficiently close to zero, so that only the t — — bands 
contribute to the integral, and 

VF± = ± \/[( & ± 6 ^) 2 ~ b cl ] i b c2 ~( b± e F ) 2 ]» ( 61 ) 

are the two Fermi velocities, corresponding to two pairs of solutions of the equation 
|b — A(k z )\ = e F , which arise from the Fermi level crossing the s = +,t = — band on 
the two sides of each Weyl node along the z-axis in momentum space. The two Fermi 
velocities are nearly equal when the band dispersion near the nodes is almost perfectly 
linear, but start to differ significantly when deviations from linearity become noticeable. 
Explicitly, as long as b c \ b ± e F b c2 , both Fermi velocities are independent of the 
Fermi energy and thus /de F vanishes. A useful way to think about this is in terms 
of an approximate chiral symmetry , which emerges in this system in the limit of small 
e F , and which manifests in an almost perfectly linear band dispersion near the Weyl 
nodes. 

By a nearly identical calculation it is easy to show that da xy /de F also vanishes 
when e F is sufficiently close to zero. This is the origin of the quasi-plateau behavior 
in Fig. [2] This result implies that the intrinsic anomalous Hall conductivity is equal 
to its thermodynamic equilibrium part, cW, not just when the Fermi energy coincides 
with the Weyl nodes, but even away from them as long as the band dispersion may be 
assumed to be linear [39]. This property will be discussed in more detail in Section |4j 

When attempting to understand the physical meaning of the two contributions to 
the anomalous Hall conductivity, a xy and a xy , it may be tempting to say that, while a : xy 
is clearly associated with states on the Fermi surface, a xy might perhaps be associated 
with the chiral Fermi arc edge states. Unfortunately, such an interpretation is clearly 
incorrect, at least in the context of the present model, since, as can be seen in Fig. 3| a T xy 
may be nonzero even when Weyl nodes and thus the Fermi arc edge states are absent. 
More work is thus needed to fully understand the relation between a xy and the edge 
states. 
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In this subsection we consider the second part of the topological response in Weyl 
semimetals, namely the CME. This arises due to the presence of an energy difference 
between the Weyl nodes, as in Eq. 34 To induce such an energy difference in our 


microscopic model, we need to find an operator that acts as an “axial chemical potential” 
term, shifting the nodes with different chirality in opposite directions in energy. The 
expression for this operator may be found based on symmetry considerations. We define 
an “axial charge density” operator, h a , as a local operator, which is odd under inversion 
and z —>■ — z reflections, even under time reversal, but odd under time reversal combined 
with rotation of the spin quantization axis by n around either x or y axis. This leads 
to the following expression for the axial charge density operator 


h n = r y a 


(62) 


Adding a term —yU a h a , where fj, a is the “axial chemical potential”, to the multilayer 
Hamiltonian, introduces an energy difference between the Weyl nodes of magnitude 

A 2 jX a V F 

Af= A 

where 

VF = 4 - (&)((& - E), (64) 

is the ^-component of the Fermi velocity at the location of the Weyl nodes. 

Similarly to the previous subsection, we now couple the electrons to electromagnetic 
field and the axial chemical potential and integrate out the electron variables to obtain an 
induced action for the electromagnetic field. We will assume, without loss of generality, 
that the electromagnetic field consists of a magnetic held in the ^-direction, and a 
vector potential A z , whose time derivative gives the ^-component of the electric held 
E z = —dtA z . We will allow for a time and ^-coordinate dependence of the vector 
potential A z and of the axial chemical potential /i a , but assume that the magnetic held 
is time-independent and uniform. 

For this calculation it is convenient to use the Landau level (LL) basis of the 
multilayer placed in a uniform external magnetic held along the growth direction, i.e. 
the basis of the eigenstates of the following Hamiltonian 

'H(k z ) = v f t z (z x <j) ■ (—zV + eA) + A (k z ) + bcr z , (65) 



Adopting Landau gauge for the vector potential A = xBy , we obtain the following 
expressions LL eigenstates, which have the form, typical for LLs in Dirac systems 

I n,a,k y ,k z ) = [z^ r (k z )\n - 1, k y , k z , t, r) + z^ T (k z )\n, k y , k z , f, r)] .(66) 


Here 


(r| n,k y ,k z ,a,r) 



( r )l^ r )> 


(67) 
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4> n k y (r) are the Landau-gauge orbital wavefunctions, and a, r are the spin and pseudospin 
indices respectively. Finally, the four-component eigenvector | z^(k z )) may be written 
as a tensor product of the two-component spin and pseudospin eigenvectors, i.e. 
| z“{k z )) = \v“{k z )) <g> | u a (k z )), where 


Vn\k z )) 

u^h)) 




m t (k z ) 

e nt(k z ) ’ 



A +A S + A D e~ ik * d \ 

V ’ A (k z ) ) 


and the eigenstate energies are given by 


e n st(k z ) = sjlu^ri + niiil:.) = se nt (k z 


( 68 ) 


(69) 


where ujb = v F /Eb is the Dirac cyclotron frequency and Eb = 1/VeB is the magnetic 
length. 

As in all Dirac systems, the lowest n = 0 LL is special and needs to be considered 
separately. The s quantum number is absent in this case and taking B > 0 for 
concreteness, we have e nt (k z ) = —m t (k z ), and \vq( k z )) = ( 0 , 1 ). 

The topological term, of interest to us, is proportional to the product of p a and A z . 
Integrating out the electron variables and leaving only this term in the imaginary time 
action, we obtain 


S = B ^ 1% itt)A z (q, itt)n a (-q , -iff) 

q,ifl 

where the response function If (< 7 , ifl) is given by: 

e n F [£ 

na' (k z )\ -n F [£ 

na (k z + q)] 


If (< 7 , iVt) = 


27 rL: 


E Uf, 

i 


n,k z 


“ 1 “ £ s na'{k , z) ^naij^z Ql) 


X , I Jz (K) l<L,) IVI z\). 


(70) 


(71) 


Here n F is the Fermi-Dirac distribution function, E, na (k z ) = £na(k z ) — Ff, and the 
magnetic field B in Eq. 70 arises from the Landau level orbital degeneracy. We have 
also ignored the g-dependence of the matrix elements in Eq. 71 which is not important 
for small q. 

At this point we will specialize to the case of an undoped Weyl semimetal, i.e. set 
e F —— 0. Then it is clear from Eq. 71 that for Landau levels with n > 1, only terms 
with s ^ s' contribute due to the difference of Fermi factors in the numerator. We are 
interested ultimately in the zero frequency and zero wavevector limit of the response 
function n(g,if2). As already seen in the previous section in the context of AHE, the 
value of n(0, 0) depends on the order in which the zero frequency and zero wavevector 
limits are taken. Below we will consider both possibilities separately and discuss their 
physical meaning. 

Before we proceed with an explicit evaluation of the q — s- 0 and O —> 0 limit, let us 
note an important property of the response function n(g, D). If we take into account 
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the following symmetry properties of the matrix elements in Eq. 71 


( V +t \v~?' 

\ u nk z I u nk z 


•i +t 
- nk , 


iv +t ' \v~ f 


), 


\cr z \v„ 


,) \ u nk z 

J nk z ) = ( V nkJ < 7 *'\ V nk z ) i O 7 ^) 

it is easy to see that when the limit 0 ^ 0 is taken, independently of the value of q, the 


n > 1 Landau levels in fact do not contribute at all, mutually cancelling due to Eq. [72] It 
follows that II(g,ffl) at small 11 is determined completely by the contribution of the two 
n — 0 Landau levels, whose energy eigenvalues and the corresponding eigenvectors are 
independent of the magnetic field. This means that in the small 11 limit II(g, 111) becomes 
independent of the magnetic field and the effective action in Eq.[70]then depends linearly 
on both p a and B, independently of the magnitude of B. 

Let us now proceed to explicitly evaluate 11(0,0), which determines the low- 
frequency and long-wavelength response of our system. Let us first look at the 
situation when we send q to zero before sending H to zero. Explicitly evaluating the 
matrix elements in Eq. [71] and the momentum integrals we obtain the following simple 
expression for II(0, 0) 

c 2 2v f 


n(o,o) = 


(73) 


47t 2 A sd 

Thus, after Wick’s rotation r —)■ it, Ae —)■ —iAe, we finally obtain the following result 
for the electromagnetic field action 

e 2 Ae 2vi 


S = 


An 2 A gd 


B / d 3 rdt A z (r,t)p a (r,t), 


(74) 


which, taking into account Eq. 63 has precisely the form of Eq. 40 Functional derivative 


of Eq. 74 with respect to A z gives the current that flows in response to magnetic field 
and axial chemical potential, i.e. the CME 
e 2 Ae 


Jz 


47T 2 


B. 


(75) 


The physical interpretation of Eq. [75] requires some care. The issue is again the 
order of limits when calculating 11(0,0), which always arises when calculating response 
functions in gapless systems. Eqs. |74|75 was obtained by sending q —> 0 before hi —> 0. 
Let us now see what happens when fl is sent to zero before taking the limit q —)• 0. 


In this case, in addition to the contribution to 11(0,0), given by Eq. 73, which arises 
due to transitions between the t — + and t = — lowest (n = 0) Landau levels, there is 
an extra contribution due to the intra-Landau-level processes within the t = — Landau 
level, which crosses the Fermi energy at the location of the Weyl nodes. This extra 
contribution is given by 

dn F (e) 

e=—m-(k z ) 


ri(o,o) = 


2i tL- 


E 


de 


(z 0kz I Jz(k z ) | z okz ) {z okz I tV I z o kz ) > ( 76 ) 


and is easily shown to be equal to Eq. [73] in magnitude, but opposite in sign, which 
means that in this case 11(0,0) vanishes. Thus, the final result for 11(0,0) depends on 
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the order in which the q —» 0 and 0 —0 limits are taken 

c 2 2v f 


lim lim II (n, O) = 

f2->0q->0 

lim lim n(o, Q) = 0. 

q->0 n-s>0 


47 t 2 A sd' 


(77) 


What is the physical meaning of these two distinct orders of limits in calculating 11(0, 0)? 
This is again identical to the AHE case, discussed in the previous section. When q is 
taken to zero first, one is calculating the low-frequency limit of response to a time- 
dependent external field, in our case magnetic field along the ^-direction. This response 
is hnite and represents CME, described by Eq. 75 If one takes hi to zero first, however, 
one is calculating a thermodynamic property, in our case change of the ground state 
energy of the system in the presence of an additional static vector potential in the 
^-direction 

1 dE(A z ) , . 

r = = 0. (78) 

This could be nonzero in, for example, a current-carrying state of a superconductor, 
which possesses phase rigidity, but vanishes identically in our case 


4. Diffusive transport in Weyl metals 

In the previous section we have focused on the basic properties of electromagnetic 
response of clean Weyl semimetals and metals (doped Weyl semimetals), neglecting 
the effect of impurity scattering. It is, however, always present, and while probably 
not of significant importance in an undoped Weyl semimetal with density of states at 
the Fermi energy vanishing as it may be expected to play an important role in a 

Weyl metal with e F significantly different from zero. If nothing else, impurity scattering 
is crucial for establishing a steady state under applied electric field, and thus must be 
included in any serious discussion of low-frequency transport phenomena. In this section 
we will thus generalize the theory of electromagnetic response, presented in Section [3j 
to the case of diffusive transport in dirty Weyl metals. As will be demonstrated below, 
this generalization leads to two important new results. 

First, we will demonstrate that the anomalous Hall conductivity of a Weyl metal 
is purely intrinsic and universal, lacking both the so-called extrinsic contribution 
due to the impurity scattering and also the part of the intrinsic electronic-structure 
contribution, coming from incompletely filled bands. Only the purely intrinsic semi- 
quantized contribution, a J J y , arising from completely filled bands and proportional to 
the separation between the Weyl nodes, survives in the Weyl metal, as long as the 
Fermi energy is close enough to the Weyl nodes, such that the band dispersion at the 
Fermi energy may be assumed to be linear. 

Second, we show that generalization of the chiral magnetic effect to the case of 
diffusive transport leads to a novel weak-field magnetoresistance effect, negative and 
quadratic in the magnetic field. We argue that this effect may be regarded as a universal 
smoking-gun transport characteristic of all Weyl metals. 
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We start again from the heterostructure model of Weyl semimetal described in detail 
in the previous sections. We would like to evaluate the anomalous Hall conductivity 
of this model ferromagnetic Weyl metal, in the presence of impurity potential V (r) = 
HoEa^r — r a ), which we will assume for simplicity to be gaussian, with only second 
order correlators present: (H(r)H(r / )) = 7 2 <5(r — r'), where y 2 = nfVf and rq is the 
impurity density. The higher-order correlators, which are known to be important for 
AHE in principle, as they give rise to the so-called skew-scattering contribution to the 
Hall conductivity, do not in fact affect our results in a significant way. We will also 
assume that the impurity potential is diagonal in both the pseudospin r and the t — ± 
index, which labels the eigenstates of the A (k z ) operator Eq. 19 Again, this assumption 
is used only for computational simplicity and does not affect the essence of our results. 
To find the anomalous Hall conductivity, we will use the method of Section 322 which 
we find to be the most convenient one for our purposes, as it allows to most clearly 
separate physically distinct contributions to the AHE. Namely, we imagine coupling 
electromagnetic field to the electrons and integrating the electron degrees of freedom 
out to obtain an effective action for the electromagnetic field, which, at quadratic order, 
describes the linear response of the system. The part of this action we are interested 
in has the appearance of a Chern-Simons term, which, adopting the Coulomb gauge for 
the electromagnetic vector potential V ■ A = 0, is given by 


q 


z°a/3 n (q, in)A 0 (-q, -in)q a Ap( q, iQ), (79) 

where q a = q a fq and summation over repeated indices is implicit. The ^-direction in 
Eq. 79 is again picked out by the magnetization b. As we will be interested specifically 
in the zero frequency and zero wavevector limits of the response function n(q, iff), we 
will assume henceforth that q = qx, which does not lead to any loss of generality due 
to full rotational symmetry in the xy- plane. The anomalous Hall conductivity is given 
by the zero frequency and zero wave vector limit of the response function n(q, ifl) as 


a 


xv = lim lim -n(q,ifi). 

y g—>0 q 


(80) 


A significant advantage of Eq. 80 compared to the more standard Kubo formula for the 
anomalous Hall conductivity, which relates it to the current-current correlation function, 
is that Eq. [80] ties the Hall conductivity to the response of a conserved quantity, i.e. the 
particle density. This means, in particular, that the response function n(q, iff) must 
satisfy exact Ward identities, which follow from charge conservation, providing a useful 
correctness check on the results. 

The impurity average of the response function n(q, ifl) may be evaluated by the 
standard methods of diagrammatic perturbation theory. Due to our assumption that the 
impurity potential is diagonal in the band index t , we can do this calculation separately 
for each pair of bands, labeled by t, and then simply sum the individual contributions. 
We will thus omit the t index in what follows, until we come to the final results. The 
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retarded impurity averaged one-particle Green’s functions have the following general 
form 


%*( k,e) = 


^ktJi ^ktJ2 


(81) 


£ - £t + 

Here s = ±, as before, labels the positive and negative energy pairs of bands (the sum 
over s is made implicit above) 

£k = se k - hF = s^Jv 2 F (kl + k 2 ) + m 2 (k z ) - e F , (82) 

are the band energies, counted from the Fermi energy e Fl and 


l^k) = —p 


1 + s 


m(k z 

Ck 


-, —ise llp \ 1 — s 


m{k z 

Ck 


(83) 


is the corresponding eigenvector. In what follows we will assume, for concreteness, that 
e F > 0, i.e. the Weyl metal is electron-doped. The impurity scattering rates l/r± are 
given, in the Born approximation, by 


1 


1 


T s (k z ) T 

where 1/r = 7ry 2 g(e F ) and 


1 + s 


m(k z )(m) 


9^f) = 


d 3 k e F 

<He k - € f) = 


t /d 


(2vr 


4n 2 v 2 F J_ v/d 


dk z O(e F - \m(k z )\), 


(84) 


(85) 


is the density of states at Fermi energy. We have also defined the average of m(k z ) over 
the Fermi surface as 


(m) = 


d 3 k 


■m(k z )6(e k - e F ). 


( 86 ) 


g(e F ) J (2 tt) 3 

The impurity averaged response function, analytically continued to real frequency as 
n(q, Q) = n(q, id —>■ kl+irj), is given, in the self-consistent non-crossing approximation, 
by the sum of ladder diagrams, as illustrated in Fig. [4] 

n(q,H) = n / (q,H) + n // (q,H), 


where 


and 


r c 

n / (q, fi) = 2e 2 v F Q / 

J —a 


de dn F (e) 
2iii de 


de 


e — irj, e + Q + irj), 


Yl n (q, Q) = Aie 2 v F / - ni?(e)ImP 0 x (q, e + irj, e + H + irj). 

J-oc 2 m 

The 4x4 matrix P, whose Ox component we are interested in, is given by 

P(q, ^ + iv) = T 2 d RA {<h fi)^(q, ft), 

P(q, e + ir 7 , e + H + irj) = I RR (e, q, Q), 

where 

V = (1 - I RA )-\ 


(87) 

( 88 ) 

(89) 


(90) 

(91) 
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n(q ,fi) 




Figure 4. Graphical representation of the impurity-averaged response function 

n(q,fi). 


is the diffusion propagator and 


tRA, 
1 ap 1 


tRR, 

1 afi ' 




« r /3 

CT2CT1'0-30-4 


e,q,fi) = -t‘ 


r P 

O-2O-1 ' O-3O-4 


(2vr) 3 CT1<T3 

d 3 k 

(2vr) 3 ffl<T3 


(k + q,O)G^ 2 (k,0), 

R (k + q, e + fi)G^ 4(T2 (k, e). 


(92) 


The physical meaning of the two distinct contributions to the response function 
n y,// (q, ff) is clear from Eqs. 88, 89 II / (q, f2) describes the non-equilibrium part of 


the response that happens at the Fermi surface, as clear from the appearance of the 
derivative of the Fermi-Dirac distribution function in Eq. [88] This response is diffusive 
when Or <C 1 and ballistic in the opposite limit, and is generally affected signihcanty 
by the impurity scattering. We will discuss this in more detail below. In contrast, 
II // (q, hi) is an equilibrium, nondissipative contribution to the overall response, to which 
all states below the Fermi energy contribute [48], and which is essentially unaffected by 
the impurity scattering, as will be seen below. 

We will start by evaluating the nonequilibrium part of the response function, 
II / (q, ff). Computation of the matrix elements q, 12) is easily done in the standard 
way, assuming f f r> 1. One obtains 


n 7 (q, Q) = ie 2 v F Q.T 9 {e F ) [I^ A V 0x + I^ A V XX 


+ 2 , 


RA V 1 

Oz ^ZX J) 


(93) 


where we have taken into account that V yx = 0 by symmetry. The relevant matrix 
elements of the diffusion propagator can be found analytically to first order in q. One 
obtains 


n J (q,0) 


ie 2 VF^lTg(e F ) 


tRA(i _ tRA\ i tRAtRA 

1 0x V 1 J zz J ' 1 0z 1 zx 

r(i - ISA) 


(94) 
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Figure 5. (a) Plot of the band edges along the ^-direction in momentum space for 
the two bands that touch at the Weyl nodes, using Ajj/A s = —0.9. (b) Field lines 
of the Berry curvature in the k y = 0 plane, for the same band structure as in (a). 
Corresponding Fermi surface section is shown by the two contours, enclosing the Weyl 
nodes. 


where 


r(q,!2) = (1 - 0(1 - o - or. 


(95) 


is the determinant of the Oz block of the diffusion propagator, which corresponds to the 
diffusion of the charge density, a conserved quantity (this block decouples from the rest 
of the diffuson when q —>■ 0). This means, in particular, that T must satisfy an exact 
Ward identity T(0,0) = 0, which is easily checked to be true. Explicitly, the relevant 
matrix elements of I RA (q, f2) to first order in q are given by 

r+/r \ RA _ iv F q /m r+/r 

1 0x ~ 


tRA _ 

-'oo — 


1 — iflr-L 


2 e F \ep 1 — ifbq 
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'0 2 
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— 1 z0 — 
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Op 1 
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jRA _ d 


iv F q 
4e f 
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1 + “ 2 “ 


2X t+/t 


1 — iQr .l 
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t + /t 


€p 1 — 


m 
— ~7T 


r+/r 
1 — iflr-L 


(96) 


where the average over the Fermi surface is defined in the same way as in Eq. 
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Expanding T(0, 12) to first order in 12 and taking the limit 12 0 at fixed r, which 

corresponds to the diffusive limit, we obtain 


0 ) — 


iqe 2 v 2 F g(e F ) / mr+\ 

——\W F[ro1 ’ 


2e, 


where 


F[m\ = 


1 + 5 


1-^) ^ 
F 


ar(o,i2) 


0 = 0 . 


-1 


(97) 


(98) 


The explicit form of the functional F[m] is in fact not that important for our purposes, 


except for the evenness property, easily seen from Eq. 98 F[m] = F[—m\. As a 
consequence, II^ is an odd functional of m, which will play an important role below. 
It is important to note that the charge conservation, whose mathematical consequence 
is the presence of the diffusion pole in II i (q, 12), is crucial in obtaining a nonzero result 
in the diffusive limit in Eq. [9T| The analogous quantity in the calculation of the spin 
Hall conductivity, for example, would vanish in the diffusive limit [521135. ;5T[[55J- 

It is also of interest to examine the ballistic limit of n 7 , which corresponds to the 
case of a clean Weyl metal. In this case we send both 12 and 1/r to zero in such a way 
that 12 t —> oo. In this case we obtain 


n L(q>0) = lira l lim n J (q,12) = lim lim ie 2 v F nrg{e F )I^ 

O-s-0 l/r—s-0 !l->0 l/r->0 

iqe 2 v F g(e F ) / m ' 

2 e F \ e F , 

which agrees with the clean Weyl metal result Eq. [56} A nice feature of the present 


(99) 


is 


calculation, compared with the calculation in a clean Weyl metal in Section 322 
that here it is particularly clear that n i (q, 12) is associated with states on the Fermi 
surface. As seen from Eqs. 


97 


and 


99 


the difference between and n£ aZ is only 


dif 


quantitative. In the AHE literature, this difference is said to arise from the so-called 
side-jump processes [52, SUEZ, SUES)- The final step of the calculation is to evaluate 
the equilibrium part of the response function, n y/ (q, 12). In the limit e F r 1 one finds 
that this part of the response function is unaffected by the impurity scattering and is 
given by 

d3 k , v, . .. a . a| scM^+q-^) ~'M£k) 

/Wk+ql r Hk/" 


n /7 (q,12) = e 2 v F J 


(2vr; 


3 \'k I ^k+q 


- « +q + H 


( 100 ) 


where summation over the band indices s, s' is again implicit. Evaluating Eq. 100 in the 
small 12 and q limit gives 


n /7 (q, 0 ) = 


—iqe 
87 r 2 


t j d 


' —n/d 


dk z sign[m(k z )] {1 - ©[^f - \m(k z )\]} 


( 101 ) 


The first term in Eq. 101 arises from the completely hllcd bands, while the second is 
the contribution of the incompletely filled bands. 
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cr xy h d/e 2 



£f/As 


Figure 6. Plot of the total anomalous Hall conductivity (solid line), a xy (dashed line) 
and a xy (dotted line) versus the Fermi energy for the same system as in Fig. [ 5 ] The 
plateau-like feature in a xy correlates with the range of the Fermi energies, for which 
the Fermi surface consists of two separate sheets, each enclosing a single Weyl node. 


We can now finally evaluate the anomalous Hall conductivity. We will focus on the 
diffusive limit results, as ballistic limit has already been discussed in Section |3j At this 
point we will also explicitly include the contribution of both t = ± pairs of bands, which 
simply amounts to restoring the index t in m t , and summing over t. Using Eqs. 80|97 


and 101 and remembering that Ao —>- iA$ upon Wick rotation to the real time, we 
obtain 


a xy 


e 2 v 2 F 


2 e, 




and 


a 


11 

xy 


= ^ ^2 I dk z sign[m t (k z )) {1 - 0[e F - | m t (k z )\)} 


r/d 


( 102 ) 


(103) 


' —ir/d 


Since m + (k z ) is positive throughout the Erst BZ, while m_{k z ) changes sign at the 
Weyl nodes, the first term in Eq. 103 which comes from completely filled bands, gives 
a universal (almost) quantized contribution 

e 2 JC 


a 


quant _ 

xy 


4(/r 2 ’ 


(104) 


where /C is the distance between the Weyl nodes, which is the same as the clean well 
metal result, Eq. 57 This equation also describes the cases when the Weyl nodes are 
absent, in which case a < ^ j ant is truly quantized since K, — 0, G, where G = 2ir/d is a 
reciprocal lattice vector. 

The most important new result of this subsection comes from examining the 
remaining, non-quantized parts of cr xy . Suppose we have a situation when the Weyl 
nodes are present and e F , while not zero, is not too far from it, as shown in Fig. [5] 
Recall that at the location of the Weyl nodes m-[k z ) = b — A {k z ) changes sign. This 

implies that, as long as /C 4^- 6f, where k ?0 is the location of a given Weyl 

z k z =k z0 

node, the average of any odd function of m_(fc 2 ) over the Fermi surface will vanish. 
This means that in such a situation, which we call Weyl metal , all contributions to 
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the anomalous Hall conductivity, associated with incompletely filled bands, will vanish 
and a xy attains a universal value, characteristic of Weyl semimetal a xy = cr^ ant , where 


chxva X y auuumu c* umv oium v mxmv,, vnui nccoiiuum wx uminuiouni kj X y ^ X y ) vv nmo 

(T X y ant is given by Eq. 104 This may be seen explicitly in Fig. [6j Note that the linear 


dispersion sufficiently close to Weyl nodes is a topological property, in the sense that 
it follows directly and exclusively from the existence of a nonzero topological charge by 
the so-called Atiyah-Bott-Shapiro construction mi. 

To understand this result physically, recall that the Weyl nodes are monopole 
sources of the Berry curvature ilk- In a clean metal, the anomalous Hall conductivity 
is given by the integral of the z-component of the Berry curvature over all occupied 

6 ( when the Fermi surface 


a 


xy 


states a xy = e 2 f ^y 3 ?i F( e k)^k- However, as clear from Fig. 
breaks up into disconnected sheets, enclosing individual nodes, the contribution of the 
states, enclosed by the Fermi surface, to this integral will always be very small, vanishing 
exactly in the limit when the band dispersion away from the nodes may be taken to 
be exactly linear. A useful analogy here is with the electric field of a dipole. A pair of 
Weyl nodes is like a dipole of two topological charges. Its field has a well-defined and 
nonzero on average ^-component at large distances from the dipole. At short distances, 
however, the field is that of individual charges, which winds around the location of each 
charge and thus any particular component of it averages to zero. 

We have thus demonstrated that the AHE in Weyl metals has a purely intrinsic 
origin and can be associated entirely with the Weyl nodes, just as in the case of a 
Weyl semimetal, when the Fermi energy coincides with the nodes and the Fermi surface 
is absent [60] . This is in contrast to an ordinary ferromagnetic metal, in which the 
anomalous Hall conductivity always has both a significant Fermi surface contribution 
and an extrinsic contribution. This property of magnetic Weyl metals may be thought 


of as being a consequence of emergent chiral symmetry, as discussed in Section 3.2 


C 2. Magnetoresistance in Weyl metals 

In this section we will discuss measurable consequences of the Chiral Magnetic Effect, 


discussed in Section 343 in the diffusive transport regime. As we will demonstrate, 
CME in this regime leads to a novel type of magnetoresistance: negative and quadratic 
in the magnetic field, first discovered by Son and Spivak |62l 63] in TR-symmetric Weyl 
semimetals. As we will show below, this novel magnetoresistance is in fact a universal 
feature of all types of Weyl semimetals and may be regarded as their smoking-gun 
transport characteristic 


We start from the axial charge density operator, given by Eq. 62 We now ask the 
following question: does h a represent a conserved quantity, as it would in a low-energy 
model of Weyl semimetal? Only when n a is conserved, or nearly conserved, may we 
expect it to contribute significantly to observable phenomena, at least at long times 
on long length scales. To answer this we need to evaluate the commutator of h a with 
the Hamiltonian "H(k). As before, it is convenient at this point to apply the following 
canonical transformation to all the operators: a ± —> t z o ± , —> a z r ± . Evaluating the 
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commutator at the Weyl node locations, we now obtain 



(105) 


This means that n a may indeed be a conserved quantity in the Weyl semimetal or weakly- 
doped Weyl metal, provided A D > A s and b = \jb c ib c2 , i.e. the magnitude of the spin 


splitting is exactly the geometric mean of its lower- and upper-critical values, at which 


the transitions out of the Weyl semimetal phase occur. Otherwise, the commutator 
is nonzero and n a is not conserved. However, as will be shown below, the relevant 
relaxation time may in fact still be long, even when the above condition is not exactly 
satisfied, in which case the axial charge density is still a physically meaningful quantity. 
The near-conservation of the axial charge density n a may be viewed as a consequence of 
an emergent low-energy chiral symmetry , which is an important characteristic feature 
of Weyl semimetals. Both the quadratic negative magnetoresistance, and the fully 
intrinsic disorder-independent AHE, discussed above, are consequences of this emergent 
symmetry. 

We now want to derive hydrodynamic transport equations for both the axial charge 
density n a (r,f) and the total charge density n(r,t). As will be shown below, what is 
known as chiral anomaly will be manifest at the level of these hydrodynamic equations 
as a coupling between n a and n in the presence of an external magnetic field. This 
coupling leads to significant observable magnetotransport effects, provided the axial 
charge relaxation time, calculated below, is long enough. As we will show explicitly 
below, long axial charge relaxation time is a direct consequence of the emergent chiral 
symmetry, which exists in a Weyl metal at low energies. 

To proceed with the derivation, we add a constant uniform magnetic field in the 
z direction B = Bz and a scalar impurity potential K(r), whose precise form will be 
specified later. Adopting Landau gauge for the vector potential A = xBy, the second- 
quantized Hamiltonian of our system may be written as 





+ 


naky k z ,n'a' k' y k' z 


Here e na (k z ) are Landau-level (LL) eigenstate energies of a clean multilayer in magnetic 
field, n = 0,1, 2,... is the main LL index, k y is the Landau-gauge intra-LL orbital 
quantum number, k z is the conserved component of the crystal momentum along the 2 - 
direction, and a = (s,t) is a composite index (introduced for compactness of notation), 
consisting of s — ±, which labels the electron- (s = +) and hole- (s = —) like sets of 
Landau levels, and t — ±, which labels the two components of a Kramers doublet of 
LLs, degenerate at b = 0. The LL eigenstate energies and the corresponding eigenvectors 
have been derived in Section 13.31 

To proceed, we will make the standard assumption, which we also made in the 
previous subsection, that the impurity potential obeys Gaussian distribution, with 
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Figure 7. Graphical representation of the impurity-averaged Green’s function within 
SCBA and of the diffusion propagator. 


(V(r)V(r / )) = 7 2 <5(r — r'). To simplify calculations further we will also assume that the 
momentum transfer due to the impurity scattering is smaller than the size of the BZ, 
i.e. | k z — k' z \d <C 1. In this case (w t (£y)|'U t, (/d.)) ~ Sue he. the t quantum number may 
be assumed to be approximately preserved during the impurity scattering. We treat 
the impurity scattering in the standard self-consistent Born approximation (SCBA), 
illustrated in Fig. [7} The retarded SCBA self-energy satisfies the equation 


yR 

^nakyk 


M 


Y (|(w,a,A: y ,A:a|y|n , ,a , ,A:(,,A:')| 2 )G^ 0 , fc j ifc ,(a;), 


n' a' k'k' 

y & 


(107) 


We will assume that the Fermi energy e F is positive, i.e. the Weyl semimetal is electron- 
doped, and large enough that the impurity-scattering-induced broadening of the density 
of states is small on the scale of the Fermi energy e F pH]. We can then restrict ourselves 
to the electron-like states with s = + (we will drop the s index henceforth for brevity), 
and easily solve the SCBA equation analytically. We obtain 


ImS 


R _ 
ntk z 


l 

2 r t (k z ) 


where 1/r = n'y 2 g(ep) and 


1 

27 


' m t (k z )(m t ) 

in - 0 


(108) 
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1 

2n£ 2 B 



nt 


(109) 


is the density of states at Fermi energy. We have also introduced the Fermi-surface 
average of m t (k z ) as 


(m t ) 


1 

27r£ 2 B g{e F ) 



y ^m t (k z )S[e nt {k z ) - e F \. 

nt 


( 110 ) 
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We are interested in hydrodynamic, i.e. long-wavelength, low-frequency density 
response of our system. As is well-known [65], the relevant information is contained in 
the diffusion propagator, or diffuson T >, given by the sum of ladder impurity-averaging 
diagrams. This is evaluated in the standard manner and we obtain 


V 1 (q, O) = 1 - /(q,G), 
where / is a 16 x 16 matrix, given by 


L a 102,0304 


(q, fi) = 


7 


LxLyL z 


( 111 ) 


d 3 r£r'e-'^>G^r,r'\n)G^,y.r\0), 


( 112 ) 

where we have introduced a composite index a = (a, r) to simplify the notation. The 
impurity-averaged Green’s functions G R,A are given by 

(r, a\n, t, k y , k z )(n, t, k y , k z |r', a') 


G R J( r,m= £ 


tt-£,nt(k z ) ±i/2r t (k z 


(113) 


where £ nt (k z ) = e nt (k z ) - e F . 

In general, the evaluation of Eq. 112 is a rather complicated task, primarily due to 
the fact that the impurity scattering will mix different LLs. At this point we will thus 
specialize to the case of transport along the ^-direction only, as this is where we can 
expect chiral anomaly to be manifest. In this case the contributions of different LLs to 


Eq. 112 decouple. Setting q = qz, we obtain 


L 0102,0304 


(q,Q) = 


7 


2tt£ 2 b L 


ST' ( a il Z n(kz + q/2)){ zi n(kz + g/2)|q 3 ) 

n “_ G - £ nt (k z + q/ 2) + i/2r t (k z + q/2) 


x 


(a 4 | Z „(k z - q/2)){z t / l {k z - q/2)\a 2 ) 


(114) 


-€nt'(k z - q/2) - i/2T[{k z - q/2) 
which is much easier to evaluate. 

As mentioned above, / and T >are large 16 x 16 matrices, which contain a lot 
of information of no interest to us. We are interested only in hydrodynamic physical 
quantities, with long relaxation times. All such quantities need to be identified, if they 
are expected to be coupled to each other. One such quantity is obviously the total charge 
density n(r,t), which has an infinite relaxation time due to the exact conservation of 
particle number. Another is the axial charge density n a (r,t), which, as discussed above, 
may be almost conserved under certain conditions. On physical grounds, we expect no 
other hydrodynamic quantities to be present in our case. We are thus only interested 
in the 2 x 2 block of the matrix D -1 , which corresponds to the coupled evolution of the 
total and the axial charge densities. To separate out this block, we apply the following 
transformation to the inverse diffuson matrix 


T >- 1 

ai 6 l,a 262 


= \(a n T bl ) 






T 7 ~ ( (T ^ T 

a2Q;i xx a 1 Q!2,^3C^4 v 


/ 0 : 30 : 4 1 


(115) 


where 01 , 2 , 61,2 = 0 ,x,y,z. The components of interest to us are 01,2 = 61,2 = 0 which 
corresponds to the total charge density, oi )2 = 0 , 6 i j2 = y, which corresponds to the axial 
charge density, and the corresponding cross-terms. 









Chiral anomaly and transport in Weyl metals 


30 


We will be interested, as mentioned above, in the hydrodynamic regime, which 
corresponds to low frequencies and long wavelengths, i.e. fir <C 1 and v F qr <C 1. 
We will also assume that the magnetic field is weak, so that ojb -C e F . Finally, we 
will assume that the Fermi energy is close enough to the Weyl nodes, so that only the 
t = — states participate in transport and (m_) ~ 0, since m_( k z ) changes sign at the 
nodes [T9]. 

In accordance with the above assumptions, we expand the inverse diffusion 
propagator to leading order in fir, v F qT and ojb/cf and obtain after a straightforward 
but lengthy calculation 




-iQr + Dq 2 r 
—iqTr 


—iqTr 

—iYlr + Dq 2 r + r/ 



(116) 


Here D = v 2 f t (m 2 _)/e 2 F is the charge diffusion constant, associated with the diffusion in 
the ^-direction, Y = eB/2ir 2 g(eF) is the total charge-axial charge coupling coefficient 
and 

1^ = 1 - (uf/A s d) 2 

r a (v f /Asd) 2 r ’ 1 ; 

is the axial charge relaxation rate. Several comments are in order here. First, an 
important thing to note is that the off-diagonal matrix elements in D _1 , proportional to 
B and responsible for the total charge to axial charge coupling, come entirely from the 
contribution of the n = 0 LL. The remaining matrix elements arise from the contribution 
of all the n > 1 Landau levels and we have taken the limit B —y 0 after summing over the 
LLs, i.e. left only the leading term in the ojb/^-f expansion. The next-to-leading term 
results in a negative correction to the diffusion coefficient, proportional to B 2 , which 
corresponds to the well-known classical positive magnetoresistance. We have ignored 
this correction here, but will comment on its effects later. Second, we note that the axial 
charge relaxation rate l/r a > 0, as it should be, and vanishes when v F = A sd. It is easy 
to see that this is identical to the condition of the vanishing of the commutator of the 
axial charge operator with the Hamiltonian Eq. |105 again as it should be. Henceforth 
we will assume that this condition is nearly satisfied so that r a r. Finally, the 
situation when v F = A s d and thus 1 /r a appears to vanish, actually needs to be treated 
with some care. Namely, the condition v F = A $d may be satisfied exactly only in the 
limit e F —y 0. The Fermi velocity depends on the Fermi energy as m 

5f(£f) = 2(6 + e F ) V / ^ + tF)2 ~^H^~ ( ' , + tF>2 l- (U8) 

When b = and thus h^(0) = A s d, the Fermi energy dependence of v F needs to 

be taken into account. Expanding to leading non vanishing order in e F we obtain in 
this case 

i 4 

T a A |T ’ 


( 119 ) 
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i.e. 1/t q is in fact always finite, but may be very small. We can estimate the minimal 
value of the axial charge relaxation rate by setting ep ~ 1/r in Eq. 119, which gives 
(■ T/T a )min « 1/(A 5 t) 2 . 

We may now write down the coupled diffusion equations for the total and axial 
charge densities, which correspond to the propagator Eq. 116 These equations read 

dn _ d 2 n dn a 
dt ~ dz 2 + dz’ 


dn a _ D ^ n a 


na r dn 
T n dz' 


dt " dz 2 r a ‘ dz' ^ 

Manifestation of chiral anomaly in these equations is the coupling between the total and 
the axial charge densities, proportional to the applied magnetic field. Since the total 
particle number is conserved, the right-hand side of the first of Eqs. 120 must be equal 
to minus the divergence of the total particle current. Then we obtain the following 
expression for the density of the charge current in the ^-direction 
cr 0 dp e 2 B 

J e dz 2n 2 /i “’ ^ ^ 

where cto = e 2 g(ep)D is the zero-field diagonal charge conductivity, pi and pi a are the total 
and axial electrochemical potentials and we have used 5n = g(ep)5p, 5n a = g(ep)5p a . 
The last relation is valid when vp/Apd is close to unity, as seen from Eq. 63 Thus chiral 


anomaly manifests in an extra contribution to the charge current density, proportional 
to the magnetic held and the axial electrochemical potential. This is known as chiral 
magnetic effect (CME) in the literature [39]. It is important to realize that the second 


term in Eq. 121 does not by any means imply that an equilibrium current may be 
driven by an applied magnetic held, despite appearances. The axial chemical potential 
/i a , appearing in Eq. 121 is a purely nonequilibrium quantity. If an equilibrium 


energy difference, p a0 , exists between the Weyl nodes due to explicitly broken inversion 
symmetry [66], then it is the difference p a ~hao that enters in Eq. 121 We have explicitly 
considered an inversion-symmetric Weyl metal here, in which case p a o = 0. 

To find measurable consequences of the CME contribution to the charge current, 
we consider a steady-state situation, with a hxed current density j flowing through the 
sample in the z-clirection. We want to find the corresponding electrochemical potential 
drop and thus the conductivity. Assuming the current density is uniform, we obtain 
from the second of Eqs. 120| 


n. 


dn 

= rr “W’ 
dz 


( 122 ) 


which is the nonequilibrium axial charge density, induced by the current and the 
corresponding electrochemical potential gradient. Substituting this into the expression 
for the charge current density Eq. |121 we finally obtain the following result for the 
conductivity 

e 4 B 2 r a 


cr = a o 


4jt 4 g(e F ) 


(123) 
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In the limit when e F is not far from the Weyl nodes, such that the dispersion may be 
assumed to be linear, we have g(ep) = e 2 F /7r 2 v F VF, which gives 

'e 2 v 2 r B' 2 


A a = cr 


Oq = 


e VFTg 

(2,71V f ') 2 


hF 


(124) 


which agrees with the Son and Spivak result [62, |63]. Thus we see that a measurable 
consequence of CME is a positive magnetoconductivity, proportional to B 2 in the limit 
of a weak magnetic held. This of course needs to be compared with the classical negative 
magnetoconductivity, which is always present and arises from the B 2 corrections to the 
diffusion constant D, which we have neglected 

A Gel 


ob 


(U C T) 


where c o c = ev 2 F B/e F is the cyclotron frequency. This gives 


Act 


Ta/r 


(125) 


(126) 


Acr rf (e F r) 2 ' 

Thus the CME-related positive magnetoconductivity will dominate the classical negative 
magnetoconductivity, provided r a is long enough. 

We have so far ignored the Zeeman effect due to the applied magnetic held. Its 
effect is to modify the spin-splitting parameter b as b —» b + gfipB j 2. In principle, 
the dependence on b does enter into our final results through the dependence of 
the Fermi velocity v F on b. Naively, this will then generate an additional linear 
magnetoconductivity, which may be expected to dominate the quadratic one at small 
fields. However, the condition of large r a , which is the same as hp/Agd ~ 1, is equivalent 
to the condition b c \ <C b <C b c2 , in which case the dependence of v F on b becomes 
negligible. Thus, in the regime in which the positive magnetoconductivity dominates 
the negative classical one, and is thus observable, one may also expect a negligible linear 
magnetoconductivity in any type of Weyl metal. 


5. Conclusions 

In this paper we have provided an overview of transport phenomena in Weyl metals, 
which may be attributed to chiral anomaly, in particular the semi-quantized intrinsic 
Anomalous Hall Effect and the Chiral Magnetic Effect, which manifests in negative 
quadratic in the magnetic field longitudinal magnetoresistance. 

The AHE is generally present in any ferromagnetic metal. In a generic FM metal 
this is a complicated phenomenon, with many physically distinct sources contributing 
to the final observed effect. In particular, extrinsic contribution to AHE, arising from 
impurity scattering, is always at least of the same order of magnitude, and often 
significantly larger, than the more theoretically appealing and universal intrinsic one. 
In a Weyl metal, however, as we have demonstrated in this paper, AHE has a purely 
intrinsic origin. As long as the Fermi energy is not too far from the Weyl nodes, such 
that the band dispersion may be taken to be linear and the chiral symmetry is present to 
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a good approximation, the anomalous Hall conductivity of a Weyl metal turns out to be 
given exactly by the semi-quantized expression of a pure undoped Weyl semimetal, where 
it simply measures separation between the Weyl nodes in momentum space in units of 
e 2 /h. Both the extrinsic, and even the part of the intrinsic contribution, associated with 
the Fermi surface (i.e. incompletely filled bands), vanish identically. We have connected 
this property with the geometry of the Berry curvature field in the vicinity of the Weyl 
nodes, which in turn is a direct consequence of their topology, through the Atiyah-Bott- 
Shapiro construction. FM Weyl metal is thus distinguished from other FM metals by 
the fact that its AHE is of purely intrinsic origin. 

Perhaps even more important in the context of Weyl metals is the Chiral Magnetic 
Effect and the associated negative quadratic magnetoresistance. This may be expected 
to occur in any type of Weyl metal, either characterized by broken time reversal 
symmetry or broken inversion (or both). It is thus of particular importance for the 
experimental characterization of Weyl metals. The negative quadratic longitudinal 
magnetoresistance should in fact be observed even in Dirac semimetals, which have 
already been realized experimentally. The theory, presented in Section |4.2 is directly 
applicable to this case, if A g = Ap limit is taken and the spin-splitting b is identified 
with the Zeeman splitting b = g/isB/ 2. In fact, this effect appears to have been observed 
very recently in a 3D Dirac semimetal material ZrTe ,5 [57]. 

What has been presented in this paper is by no means a complete story of 
the possible observable manifestations of chiral anomaly in Weyl metals. Other 
related phenomena, such as linear high-field magnetoresistance [68], plasmon-magnon 
coupling (69], nonlocal transport [70j, anomalous density response cn, anomalous 
thermoelectric response [72] , and others [721 EH ESI ESI EH ESI ESI ESI EH ESI ESI El] , 
have been discussed in the literature. One of the possible promising directions for future 
research in this area is the interplay of chiral anomaly and superconductivity in Weyl 
metals 
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